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ABSTRACT 


Coastal  Hydrographies  Techniques  (CHT)  personnel  in 
support  of  the  Airborne  Bathymetry  System  ( ABS j  traveled  on 
two  separate  occasions  to  the  CERC  FRF  at  Duck,  N  C.  ,  to 
provide  ground  truth  data  for  the  overflights  of  the  system 
The  following  data  report  gives  a  preliminary  look  at 
the  data  and  provides  some  preliminary  analysis  and  simple 
calculations  using  the  collected  data  for  both  CHT  and  ABS 
appl ications. 
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INTRODUCTION 


The  Defense  Mapping  Agency  (DMA)  has  provided  support  for 
the  deveiopement  of  the  Airborne  Bathymetry  System  (ABS)  Five 
flights  for  fiscal  year  1988  were  planned  by  NORDA  in  conjunction 
with  DMA  and  the  Naval  Oceanographic  Office  (NAVOCEANO)  These 
test  flights  of  the  system  were  to  take  place  at  selected  areas 
chosen  for  their  water  optical  quality,  bottom  type  and 

reflectivity,  and  bathymetry  database  availability. 

A  site  selected  for  the  ABS  overflights  was  the  Army  Corps 
of  Engineers  Coastal  Engineering  Research  Center  (CERC)  Field 
Research  Facility  (FRF)  at  Duck,  North  Carolina,  This 
ocean  enqineering  research  center  provides  an  ideal  area  fcr  test 
overflights  due  not  only  to  the  water  optical  quality  variability 
but,  also,  because  the  facility  research  mandate  provides 
valuable  supporting  data  tor  unatrsianoiiig  ii.c  impact  that 
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environmental  factors  have  as  noise  sources  in  the  return  signals 
of  the  HALS  laser  sounder.  The  CERC  facility  routinely  measures, 
collects,  and  archives  wave  spectra  and  direction,  prevailing 
currents,  weather  conditions,  bottom  sediment  types,  high 
accuracy  bathymetry,  tidal  information,  and  sediment  transport 
Appendix  A  provides  ancillary  information  about  CERC  FRF  and  the 
archival  database  present 

To  support  its  investigations,  the  FRF  utilizes  a  6.1  by  661 
meter  long  p'er  that  extends  from  behind  the  dune  line  seaward  to 
about  the  ?  meter  depth  contour  At  the  terminus  of  the  pier  is 
stationed  an  air-conditioned  van  which  can  be  used  to  house 
instruments  for  experiments  conducted  near  the  end  of  the  pier 
Locations  on  the  pier  are  enumerated  in  feet  from  a  baseline 
located  landward  of  the  laboratory  building  and  normal  to  the 
pier  centerline  The  laboratory  building  includes  offices,  a 
kitchen,  library,  computer  center,  a  multipurpose  area,  and  a 
diving  locker  (1). 

In  two  separate  ABS  field  support  excursions  to  the  CERC  FRF 
at  Duck,  data  was  collected  for  the  purpose  of  providing  ground 
truth  for  the  P-3  overflights  The  first  field  study  ran  from  14 
June  1988  to  19  June  1988  inclusive  with  the  overflight  occurring 
on  18  June  1988.  The  second  field  study  ran  from  14  August  1988 
to  18  August  1988  inclusive  with  the  overflight  occurring  on  17 
August  1968. 
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FIELD  INSTRUMENTATION  AND  DATA  COLLECTED 


The  instrumentation  used  for  determining  the  ocean  optical 
properties  was  a  spectroradiometer  that  was  designed  and  built  by 
Research  Support  Instruments  (RSI|.  The  instrument  measures 
irradiance  using  a  cosine  collector  at  five  channels  of  radiant 
energy  input  Four  of  the  five  channels  are  equivalent  to  the 
channels  of  the  Coastal  Zone  Color  Scanner  (CZCS)  One  channel 
is  at  a  spectral  benchmark  that  allows  interpolation  of  the 
measurements  made  at  this  wavelength  to  other  wavelengths  (see 
Appendix  B  for  wavelength  values  and  comparison  to  the  CZCS 
wavelengths].  The  spectroradiometer  data  is  read  from  the  panel 
meters  in  the  control  box.  Upwelling  and  downwelling  irradiance, 
pitch,  roll,  and  water  temperature  are  recorded  at  distinct 
depths  increments. 

Optical  measurements  were  collected  at  three  separate 
stations  on  the  first  trip  and  two  stations  on  the  second  trip 
(with  readings  taken  twice  at  one  station).  These  data  were 
taken  along  the  pier,  seaward,  with  the  RSI  instrument.  Table  1 
provides  relevant  station  information.  The  data  collection 
procedure  involved  utilizing  a  rigid  support  frame  that  allowed  a 
block  to  be  hung  over  the  handrail  of  the  pier.  Through  the 
block  ran  the  support  line  and  electrical  cable  for  the 
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instrument  The  spectroradiometer  was  lowered  to  a  position  just 
above  the  water  surface  and  initial  readings  collected  These 
readings  consisted  of  the  up  welling  irradiance  and  the 
downwelfing  irradiance.  Once  completed,  the  instrument  would  be 
placed  below  the  surface  of  the  water  and  a  second  set  of 
readings  wouid  be  recorded  From  these  data,  the  spectral 
transmission  of  the  inte-face  could  be  estimated  Moreover, 
spectral  estimations  of  the  surface  reflectance,  w  a  t  e r  leaving 
irradiance,  and  absorption  of  the  surface  water  layer  could  be 
made  Once  the  second  set  of  readings  was  completed,  the 
Instrument  was  lowered  incrementally  until  near  the  bottom  From 
these  data,  the  diffuse  attenuation  coefficient  can  be  reckoned 
and  the  spectral  bottom  reflectivity  estimated.  Bathymetric 
information  on  the  station  was  recorded.  The  three  stations 
taken  were  at  660,  1220,  and  1900  feet  respectively  for  the  first 
data  set  The  second  data  set  comprises  data  taken  at  1220  (two 
datasets  on  the  same  day  at  different  times)  and  1900  feet  on  the 
pier. 

The  primary  interest  the  water  optical  data  has  is  the 
impact  the  multifarious  environmental  factors  might  have  on  the 
active  portion  of  the  ABS,  the  HALS  laser  sounder,  return 
signals  Additionally,  these  optical  measurements  could  be  used 
to  improve  the  bathymetry  algoritms  used  in  the  passive  subsystem 
of  the  ABS,  the  NORDA  scanner. 
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DATA  REDUCTION  AND  DISPLAY 

The  data  reduction  procedure  for  each  plot  In  the  following 
set  of  graphs  will  be  briefly  discussed. 

The  spectral  transmission  of  the  interface  was  estimated  by 
looking  at  the  downwelling  above  the  surface  versus  the 
downweiling  below  the  surface  Because  wave  action  alternately 
increases  and  decreases  the  overlying  water  layer  thickness,  the 
measurement  is  somewhat  difficult  to  make  due  to  the  rapidly 
varying  readings  obtained.  Nonetheless,  by  looking  at  the  maxima 
and  minima  readings  obtained  and  taking  an  intervening  reading, 
some  reasonable  approximation  can  be  obtained.  It  is  to  be  noted 
that  the  layer  of  intervening  water  between  these  measurements  is 
not  infinitesimally  thin.  Thus,  some  absorption  invariably 
occurs.  No  provision  for  correcting  the  transmission  for  this 
absorption  is  made  in  the  relevant  plots.  Figure  1  exhibits  the 
spectral  transmission  of  the  interface  for  the  three  pier 
stations  for  the  first  Duck  field  excursion.  Relevant  data  for 
second  Duck  field  study  was  not  recorded. 

The  reflectivity  of  the  naviface  could  be  estimated  by 
comparing  the  downwelling  and  the  upwelling  irradiance  while  the 
instrument  is  above  the  surface  Care  must  be  taken  here  since 
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the  instrument  ought  not  shadow  itself  while  in  the  process  of 
extracting  a  reading  Moreover,  as  with  the  transmission 
measurement  aforementioned,  the  collimated  component  represented 
by  the  direct  sun  will  affect  the  measured  values.  That  is,  the 
transmissivity  and  reflectivity  of  the  naviface  is  dependent  on 
the  solar  angle  Theoretically,  one  could  subtract  the 
transmissivity  from  unity  and  obtain  a  value  for  reflectivity 
simply  through  conservation  of  energy  considerations  [provided  no 
absorption  occurs)  However,  emanating  from  below  the  water 
surface  through  the  air-water  interface,  is  light  due  to  the 
water  column  This  water  leaving  irradiance  is  embedded  in  the 
surface  reflectance  flux  and,  thereby,  is  part  of  that  which  is 
measured  as  reflectivity.  Figure  2  shows  the  reflect. vity  of 
the  interface  for  the  first  field  data  set.  Figures  3  and  A  show 
the  percentage  of  water  leaving  irradiance  for  the  first  data  set 
and  for  the  third  station  of  the  second  data  set  when  compared  to 
the  overall  light  leaving  the  air-water  interface 

If  it  is  certain  that  some  absorption  took  place  in  the 
interface  region,  then  since  the  sum  of  reflectivity, 

transmissivity,  and  absorption  ought  to  equal  unity,  one  is  able 
to  estimate  the  amount  of  absorption  in  the  surface  layer  of  the 
interface.  Figure  5  shows  this  kind  of  calculation  for  the  first 
Duck  data  set 

The  irradiance  reflectances  were  measured  for  each  level 
taken  at  each  station.  Notice  at  the  top  station,  above  the 
interface,  this  amounts  to  measuring  the  reflectivity  of  the 


interface  At  the  bottom,  if  such  a  measurement  could  reasonably 
be  made,  this  data  would  allow  the  reckoning  of  the  bottom 
reflectivity  The  air-water  interface  and  the  benthic  interface 
are  actual  (real)  interfaces.  The  intervening  ratios  of  the 
upwelling  to  downwelling  are  not  of  real  interfaces,  therefore, 
they  will  be  referred  to  hereafter  as  irradiance  ratios  The 
irradiance  ratios,  treating  channel  and  depth  as  the  independent 
variable  by  turns,  are  shown  for  each  station  in  Figs  6-20 
The  upwelling  and  downwelling  diffuse  attenuation  coefficient  (or 
' k ' )  was  computed  as  a  function  of  depth  and  channel.  Figures  21 
through  44  show  these  results. 

The  bottom  reflectivity  measured  by  the  laboratory 
ref lectometer  on  the  bottom  samples  taken  at  Duck  for  the 
different  stations  is  shown  in  Tables  2  and  3  The  bottom 

reflectivity  as  calculated  by  a  simple  bulk  properties  model  and 
by  linear  regression  of  the  irradiahce  ratios  as  a  function  of 
depth  for  a  given  wavelength  are  compar  ~d  in  Figs.  45  through  50 
for  the  Duck  study.  The  bulk  properties  model  is  provided  as 
Appendix  C.  The  weakness  of  the  model  is  in  assuming  the  diffuse 
attenuation  coefficients,  or  k  values,  for  both  the  upwelling  and 
the  downwelling  light  streams  for  the  last  level  above  the 
bottom,  to  be  adequate  for  the  calculation. 
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DISCUSSION  OF  DATA 


The  discussion  that  follows  is  abbreviated  and  will  be 
expatiated  in  a  follow  up  report  In  the  following  discussion, 
the  term  'channel'  refers  to  the  specific  wavelengths  of  the 
sensor  filters  possessed  by  the  spectroradiometer  See  Appendix  B 
for  details 

In  the  discussion  above.  Fig  l  was  shown  to  exhibit  the 

transmission  of  the  interface  For  a  calm  interface,  a 

calculation  may  be  made  using  the  well  known  Fresnel  formulae 
If  we  assume  a  refractive  index  of  1.34  for  sea  water,  Figs.  51 

and  52  give  the  results  for  various  angles  of  incidence  in  air 

and  water  respectively  Notwithstanding,  a  calm  interface  is 
hardly  realistic  Cordon  (2)  has  given  the  time  averaged 
reflectances  for  various  windspeeds  from  0  to  19  m/sec.  Fig  53 
gives  the  transmittance  of  the  air-water  interface  as  a  function 
of  windspeed.  Cordon  (2)  assumes  that  the  crosswind  and  upwind 
slope  parameters  may  be  represented  by  a  single  Caussian 

distribution  of  slopes  This  is  a  simplifying  assumption 
Comparing  Fig  1  results  to  Fig  51  shows  obvious  discrepancies 
For  an  observation  angle  of  zero  degrees,  the  transmittance  ought 
to  be  a  little  less  than  98  %  However,  Fig  1  provides  smaller 
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values  Furthermore,  since  there  is  some  reflection  from  the 
naviface  for  light  coming  up  from  below,  this  gets  summed  with 
the  light  the  downwelling  sensor  on  the  RSI  instrument  sees 
transmitted  through  the  interface  This  adds  a  small  positive 
bias  to  the  measurement  of  the  transmission  through  the  air-water 
boundary.  The  plots  in  the  graph  shown  in  Fig  1  have  not  been 
corrected  for  the  bias  mentioned 

Fig  2  exhibits  the  reflectivity  for  the  first  Duck  data  set 
Again,  the  values  expected  from  Gordon  (2)  do  not  compare  well 
with  the  measured  values.  For  a  nadir  observation  angle,  the 
reflectance  of  the  naviface  ought  to  be  a  little  over  2  % ,  even 

for  windspeeds  up  to  16  m/sec 

Figs  3  and  A  show  the  water  leaving  irradiance  (water 

column  light)  for  the  three  stations  associated  with  the  first 
Duck  data  set  and  for  the  third  station  of  the  second  Duck  data 
set  An  interesting  result  here  is  the  apparently  high 
percentage  of  'embedded'  flux  or  water  leaving  irradiance  in  the 
total  light  that  leaves  the  interface  upward.  Station  two 
exhibits  an  improbable  excess  of  water  leaving  irradiance 
Three  factors  may  enter  here  and  in  other  data  results  to  cause 
these  kind  of  ernneous  readings  First,  there  can  be  simple 
reading  errors  in  recording  the  measurement  from  the  panel 
meters,  or  simply  errors  involved  in  choosing  a  nonrepresentative 
reading  when  the  readings  on  the  panel  meter  fluctuate  rapidly. 
Secondly,  there  was  no  deck  cell  coupled  with  the  instrument  to 
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allow  a  way  to  normalize  the  changing  lighting  conditions  extant 
while  the  measurements  were  being  made.  Thirdly,  the  instrument 
is  lowered  into  the  water  below  the  surface  far  enough  to  have 
the  instrument  covered  completely  even  with  the  presence  of 
waves.  The  thickness  of  the  water  layer  is  therefore  a  function 
of  the  wave  action  present.  Since  there  is  a  water  layer  to 
propagate  across,  the  readings  found  for  the  upwelhng  below  the 
interface  may  not  include  the  attenuation  near  the  surface  due  to 
suspensoids,  and  transmission  through  the  interface  However, 
the  total  light  off  the  water  surface  will  include  such  effects 
The  plots  of  the  reflectivity  seen  in  Fig  2  are  not  corrected 
for  the  water  leaving  irradiance  fraction 

Fig  5  provides  an  estimation  of  the  absorption  of  the  water 
layer  when  combined  with  the  transmission  and  reflectivity  of  the 
interface  If  we  know  the  transmission  and  reflectivity,  then 
the  absorption  can  be  found  An  interesting  result  here  is  that 
station  one  in  the  first  Duck  data  set  exhibits  negative 
absorption  Or,  in  other  terms,  light  seems  to  be  created  in  the 
water  layer  referred  to.  Beyond  simply  an  error  being  made  in 
taking  the  reading,  this  possibly  may  be  due  to  the  presence  of 
suspended  particulate*  in  the  water  layer  bracketed  by  the 
instrument  in  its  first  two  measurement  positions  at  the  station. 
Station  one  was  quits  close  to  the  beach  area  This  allowed  the 
water  to  have  a  quantity  of  resuspended  sediments  near  its 
surface  layers  These  suspensoids  allowed  for  a  good  deal  of 
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backscattered  light  tc  be  seen  r>y  the  downlooking  sensor  above 
the  interface  and,  yet,  not  picked  up  by  the  uplooking  sensor 
below  the  water  surface  layer.  This  tended  to  inflate  the 
reflectivity  measured  and  reduce  the  apparent  transmission  Fig 
54  displays  the  excess  percentage  of  irradiance  at  station  one 
over  TOO  %.  Figs.  6  through  20  provide  the  irradiance  ratios 

for  both  Duck  data  sets  as  a  function  of  both  channel  and  depth 
There  were  some  apparently  aberrant  measurements  made  in  these 
data.  For  example,  Fig  18  shows  the  irradiance  ratios  for 
station  three  of  the  second  Duck  data  set  Here,  the  irradiance 
ratios  exceed  unity  Ostensibly,  the  recorded  values  were  in 
error. 

Figs  21  through  44  provide  the  upwelling  and  downwelling 
diffuse  attenuation  coefficients  (k's]  as  a  function  of  channel 
and  depth  By  and  large,  the  upwelling  k’s  exceed  the 
downwelling  k's  in  value  Presumably,  if  the  water  mass  were 
uniform,  the  upwelling  and  downwelling  k's  would  be  identical 
Thus,  the  nonequivalence  of  these  provides  a  rough  gauge  of  the 
nonuniformity  (stratification]  of  the  water  column. 

Figs.  4$  through  50  provide  estimations  of  the  bottom 
reflectance  in  two  ways  First,  the  irradiance  ratios  as  a 
function  of  depth  are  plotted  and  fit  by  a  curve  which  allows  a 
prediction  of  the  bottom  reflectance,  for  a  particular  channel, 
for  the  given  depth  Secondly,  the  water  mass  is  modelled  and 
using  as  input  the  measured  upwelling  and  downwelling  k’s  as  well 
as  the  depth  of  the  station,  the  bottom  reflectance  is  calculated 
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for  the  channel  in  question.  Again,  some  of  the  values  obtained 
by  these  methods  are  simply  too  large.  For  instance,  Fig  49 
exhibits  the  bottom  reflectivity  estimation  for  station  two  of 
the  second  Duck  data  set  Obviously,  the  values  are  too  high 
Lyzenga  (3)  provides  some  bottom  reflectivity  values  as  a 
function  of  wavelength  for  sand,  silt,  shoal  grass,  and  turtle 
grass  from  the  Panama  City,  Florida  area.  Fig  55  shows 

Lyzenga's  results  For  a  relatively  light  colored  sand  shown, 
the  percentage  of  reflectance  never  exceeds  30  % 

Tables  2  and  3  provides  the  data  sheets  for  the  measured  values 
of  the  reflectivity  of  bottom  samples  obtained  at  Duck  for 
different  distances  along  the  pier.  Certain  distances  match 
closely  the  three  distances  at  which  stations  acquired  data  and 
the  other  distances  provide  some  insight  into  the  nature  of  the 
change  of  bottom  reflectivity  as  one  moves  seaward  from  the 
laboratory  building  The  measurements  were  made  in  the 

laboratory  using  a  Hunter  Reflectometer  Model  40D  that  uses 

either  (or  both)  a  green  and  blue  filter  to  make  ’whiteness' 
and/or  reflectance  measurements.  Some  agreement  in  the  calculated 
and  measured  values  can  be  seen  by  comparing  the  results  for 
station  1.  The  filters  used  by  the  Hunter  Reflectometer  are  wide¬ 
band.  Thus,  channels  1  and  2  will  be  covered  by  the  ’blue' 
filter  and  channel  3  covered  by  the  ’green’  filter  Fig  45 
shows  the  bottom  reflectivity  for  station  1.  Comparing  channels 
1  and  2  to  the  laboratory  measurements  in  Fig  56  for  the  700 


foot  mark,  the  agreement  is  evident  for  the  extrapolation  of  the 
irradiance  ratios  to  the  bottom  However,  the  respective  values 
of  channel  3  not  only  show  different  trends  but  have  „ o m e w h a t 
different  values  The  top-most  plot  In  Fig  56  shows  the  result 
of  the  Hunter  -  Judd  whiteness  equation.  This  gives  a  direct 
proportion  scale  from  0  to  100  for  'whiteness'.  Interestingly 
enougn,  the  darker  sediments  at  the  end  of  CERC  pier  show  a 
higher  value  of  'whiteness'  than  those  closer  to  the  shore.  It 
appears  something  is  awry  here 

The  disagreement  between  the  measured  and  calculated  va'ues 
either  by  extrapolation  or  modelling  possess  can  be  due  to  the 
implicit  assumptions  that  could  be  easily  violated  by  a  complex 
water  mass  For  example,  in  the  extrapolation  process,  it  is 
expected  that  the  irradiance  ratios  vary  smoothly  with  depth 
This  is  obviously  not  the  case  Hence,  any  extrapolation  taken 
is  done  so  at  some  risk  In  the  bulk  properties  model  alluded  to 
earlier  on,  there  is  the  clear  indication  that  the  k's,  both 
upwelling  and  downwelling,  do  not  represent  well  enough  the 
actual  attenuation  which  holds  in  the  vicinity  of  the  bottom 

In  order  to  carry  the  analysis  a  somewhat  further  and 
provide  connection  to  the  ABS  HALS  laser/receiver  subsystem, 
Appendix  D  gives  ;ome  relevant  laser  power  and  signal-to-noise 
calculations  using  the  data  found  at  station  3  for  the  Duck  data 
set  1  Appendix  E  uses  an  algorithm  due  to  Austin  and  Petzold 

(4)  to  calculate  the  k  value  for  the  upper  surface  of  a  water 
mass  using  a  ratio  of  the  water  leaving  irradiance  at  two  CZCS 
wavelengths  This  k  value  is  then  compared  to  the  measured  k 


value  (see  Figure  58)  Append-x  f  exhibits  tne  field  data  sheets 
themselves 
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FUTURE  ANALYSIS  EFFORTS 

The  data  collected  by  CHT  personnel  at  Duck  CERC  FRF,  when 
used  in  conjunction  with  the  results  from  the  overflights  of  the 
ASS  will  allow  analysis  which  will  address  the  effects  of  the 
environment  on  the  ABS  signals  Beyond  the  ABS  requirements 
analysis,  further  investigations  into  the  applicability  of  the 
CZCS  algorithm  developed  by  Austin  and  Petzold  for  coastal  waters 
would  be  very  useful. 

An  aspect  touched  on  in  the  above  has  been  the  reflectivity 
of  the  bottom  sediments  collected  at  Duck  Both  laboratory 
measurements  and  some  theoretical  computations  were  utilized  to 
obtain  values  for  the  reflectivity.  It  appears  that  very  little 
actual  work  has  been  done  on  bottom  sediment  reflectivity.  Yet, 
it  is  a  very  important  parameter  in  the  prediction  and  actual 
operation  of  many  ocean  optical  systems  in  coastal  waters. 
Further  analysis  in  this  direction  is  very  much  needed. 
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personnel  at  Duck,  N.C.,  who  mace  our  -field  worn  pleasant  and 
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The  CZCS  and  seven  other  discrete  sensors  were 
incorporated  into  the  Nimbus-7  satellite  launched  from  Cape 
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FIG  1. 

Assumptions 

1.  The  sea  is  passive. 

2.  k  is  not  necessarily  constant.  It  can  change  at  each  level. 

3.  Levels  can  be  of  different  thicknesses. 
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-from  bottom  up. 
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Cal cul at i ons 


where  i  in  the  -first  exponential  goes  -from  1  to  n,  or  the 
total  number  of  layers  present.  In  the  numerator ,  i  goes  from 
1  to  1  and  in  the  denominator  1  goes  from  1  to  1 . 
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exp 
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■  exp 
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J  di  i 
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where  i  in  the  first  eponential  goes  from  1+1  to  n  and  in  the 
second  goes  from  1  to  1 


Note  that  the  first  and  second  exponential  expressions  cover 
the  same  region.  Whereupon, 


R  :=  R  exp 

_ O  • 

k  ■  u 

1  b 

ui  i 

m 

where  we  assume  the  upwelling  and  downwelling  k's  are  equal. 


In  terms  of  E  sub-u  and  E  sub-d; 
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ul  ' 
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dl 
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R  exp 
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\  ' 

/ 

k  •  u 

b 

^ _ l 

ui  i 

i 

The  best  reckoning  of  R  sub-b  will  be  when  we  are  one  interval 
above  the  bottom;  then,  l'=l  and 
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LASER  POWER  AND  SNR  CALCULATIONS 


POWER  RECEIVED 


The  following  is  a  quick,  heuristic  derivation  of  a 
laser  received  power  equation  and  a  SNR  equation  for  the 
system  Bulk  properties  of  the  water  mass  are  used 

If  the  power  produced  by  the  laser  is  P,  then  after 
exiting  the  aircraft,  the  power  has  become  diminished  by 
passage  through  the  transmitting  optics,  which  gives 

r% 

and  when  coupled  with  losses  in  the  atmosphere,  we  have  at 
the  surface  of  the  water 


TT,1 


t  'ft 


Losses  in  the  air-water  interface  are  allowed  for  by  a  factor 
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or,  altogether, 


The  laser  beam  produces  a  spot  on  the  water  surface  and  on 
the  ocean  floor  The  spot  size  on  the  ocean  floor  will  be  a 
function  of  the  beam  divergence,  turbidity  in  the  water  and 
the  water  depth  Civen  the  beam  divergence  and  the  water 
depth,  one  can  estimate  the  bottom  spot  size  using  results 
provided  by  Duntley  (5).  Duntley  provides  a  factor  to  be 
included  in  the  calculation  that  taken  account  of  the  spot 
size  and  how  it  affects  the  irradiance  at  depth  Using  this 

factor,  w,  the  reflectivity  from  the  bottom,  and  the 
attenuation  to  the  bottom  (in  exponential  form),  we  can  write 


Pt  T;Ta^  e  W  [-  Kv/cos  <Tj 


The  reflected  laser  light  traverses  the  water  on  its 
return  path  and  exits  the  surface.  It  travels  through  the 
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atmosphere  Including  the  attenuation  factor  and 
transmission  factors,  we  have 


I f  we  ass  ume 
in  the  upper 
and  decrease 


the  bottom 
hem i sphere 
as 


is  a  Lambertian  reflector,  rho,  then 
the  upwelling  radiation  will  spread 


TT  #■  -lo  rece-wt^ 


The  receiver  will  intercept  the  beam  and  degrade  its  power 
through  its  receiving  optics.  Taken  altogether,  then,  we 
have 
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^ W A  '  COSrj= 


The  angle  phi  which  represents  the  angle  in  the  water  can  be 
transformed  into  the  incidence  angle,  theta,  in  air  via 
Snell's  Law.  If  this  is  done,  the  resultant  expression  is. 


Equation  [1]  provides  a  calculation  for  the  power  received  by 
the  receiver. 


SNR  CALCULATION 


The  primary  noise  source  is  the  sun.  The  solar 
irradiance  can  be  written 
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b5  * 


where  is  the  solar  spectral  irradiance 

of  the  sun  off  the  water  can  be  given  by 


The  reflectance 


B5  *  £>\  x  [  I  -  T^) 


As  before,  the  light  is  attenuated  in  its  passage  through  the 
atmosphere  by  both  aerosols  and  geometry.  If  we  wish  to 
state  the  power  seen  at  the  receiver  due  to  the  solar 
irradiance,  then  we  can  multiply  the  spectral  solar 
irradiance  by  the  area  being  viewed.  This  involves  the  field 
of  view  of  the  receiver,  the  distance  from  the  receiver  to 
the  area  encompassed  by  the  field  of  view,  and  a 

l 

term  as  an  areal  weighting  factor.  Altoather,  including  an 
aperture  term  for  the  receiver,  we  con  write 


130 


However,  solar  light  will  also  pass  through  the  water 
surface  Thus,  another  term  needs  to  be  included. 
Essentially,  this  follows  earlier  arguements  with  apposite 
transmission  factors,  exponential  attenuation  factors,  and 
attenuation  due  to  geometric  spreading  in  the  upper 
hemisphere.  Altogether,  the  term  can  be  written  as 
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TT 


[3] 


The  total  solar  noise  becomes  the  sum  of  [2]  and  [3] 


Another  source  of  noise  is  that  due  to  the 
photomultiplier  (PMT).  The  noise  equivalent  power  (NEP)  for 


a  typical  PMT  can  be  written 


In  the  above,  N  is  the  NEP  in  watts,  q  is  the  charge  on  the 
electron,  I  is  the  anode  dark  current  in  amperes,  C  is  the 
amplification  ratio  for  the  current,  f  is  the  electronic 
bandwidth  in  Hz,  and  S  is  the  anode  radiant  sensitivity  in 
amps/wat  t  at  a  given  channel  . 


The  cal 

c  u  1 

a  t  1 

ion  of  the  signal  to 

noise  ratio 

[SNR)  is 

best 

per  f  ormed 

i  n 

photons.  One  may  t 

ransform  wat 

t  s  into 

phot 

o  ns / se  c 

by 

us 

ing  the  relation 

X 

he 

Then,  knowing  the  time  interval,  t,  during  which  our 
measurement  is  made,  and,  if  we  know  the  quantum  efficiency 
of  the  device,  then  we  can  translate  power  into  photons. 

For  example,  the  number  of  photons  available  from  the  laser 
pulse  will  be 


C5]  . 


/ 


The  noise  present,  as  previously  considered,  consists  of 
three  parts  In  terms  of  photons,  these  may  be  written  as  in 
the  following 

j 

J  IP  •  IP  ■ 

I 

These  components  can  be  rati oed  appropriately  [assuming 
Poisson  statistics)  to  give 


_ ;v _ 

CSr  +  S"  +  Y 

[6]. 


For  the  case  of  the  data  provided  by  station  3  for  the 
first  Duck  data  set,  the  following  table  gives  the  power 
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received  and  the  signal  to  noise  ratio  for  the  assumed 
parameters  In  this,  the  noise  assumed  for  the  PMT  is 
15  watts  and  the  quantum  efficiency  is  15  %  Obvious! 
other  calculations  can  be  made  if  one  wishes  to  specul 
different  system  or  water  column  parameters. 


SYSTEM  PARAMETERS 

POWEROUT  =  75  mJ/PULSE 
APERT  URE  =  6  INCHES 
MEASURE  T I  ME -2 . 5E -  09 
RECEIVER  F OV  =  008  RADS 
TRANSMISSION  OPTICS*  .5 
RECEIVER  OPT ICS*  . 5 
ATMOS PH . TRANSM 1 SS 1  ON  =  .98 

FILTER  BANDWIDTH*  . 0  1 
SCAN  ANCLE* 15  DECREES 
SPECTRAL  SOLAR 
ALT  1 TUDE  =  500m 


WATER  COLUMN  PARAMETERS 
(STATION  3  DATA) 

DIFFUSE  ATTENUATION  COEFF 
BEAM  SPREAD  FACTOR*  .16 
INTERFACE  TRANSMI SSION* . 6 
BOTTOM  REFLECTIVITY*  .  12 


micrometers 


IRRAD  =8  DOE-  02 


The  diffuse  attenuation  coefficient  used  will  be 
average  value  over  the  water  column  for  the  channel  c 


system 

1  9  5  E  - 

t  e  on 


=  577 


an 

I  oses  t 
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to  the  doubled  Nd.YAG  green  line  (532  nm )  Similarly,  the 
bottom  reflection  and  interface  transmission  will  be  taken 
from  the  data  of  the  channel  that  most  closely  approximates 
the  532  nm  line. 

POWER  RETURNED  TO  RECEIVER  SIGNAL  TO  NOISE 

7 . 595E-09  Watts  .  2564 

Obviously,  we  would  not  be  aule  to  see  the  return 
signal.  This  is  in  line  with  the  'KD  '  limits  of  the  system 
The  KD  value  for  this  set  of  parameters  is  3.17  Since  the 
limits  of  the  system  are  KD=3.00,  the  system  should  not  be 
expected  to  give  a  readable  signal  back 

Interestingly  enough,  if  in  the  above  table  of  water 
column  values  we  change  only  the  bottom  reflection 
coefficient  to  a  higher  value,  say  .3,  the  SNR  becomes  .636 
which  is  approaching  the  unity  b r eak - e ve n - po i n t  for  the  same 
KD  value  as  the  first  case. 

As  a  last  situation,  suppose  all  else  is  the  same  as 
in  the  tabulated  list  but  for  the  k  value,  which  becomes  4 
The  KD  value  of  the  system  is  now  2.2  and  the  SNR  is  3.6  --  a 
measureable  returning  irradiance. 

Fig.  57  gives  a  plot  of  k  vs  d  for  a  kd=3.0 
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APPENDIX  E 


CZCS  ALCORITM  CALCULATIONS 


« 


A  u  s  t  i  n 

and 

p  e  t  z  o I d  (4 

1  developed  an  algor 

i  t  hm 

f  o  r 

t  h  e 

c  a i c  u i a  t ion 

o  f 

the  diffuse 

attenuation  coeffic 

ient 

f  o  r 

t  wo 

the  CZCS  wavelengths  (-490  and  520  nm)  by  rat  icing  certain 
bands  of  CZCS  and  comparing  these  to  field  data  and  obtaining 
a  regression  relationship  between  the  data  and  the  ratio 


The  relevant  algorithm  for  the  490  nm  k-value  can  be 
written  as 


fc(4o,o)-  0.0SY 3 


L(44 3)  v-mh. 

l~G wj 


+ 


[1  ]  • 
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The  germane  algorithm  for  the  520  nm  band  is 


K(s2®)=  d>.Q>U3 


The  graph  in  Fig  58  gives  the  calculated  values  of  k 
for  both  490  nm  and  520  nm  to  compare  with  the  measured  k  for 
the  stated  stations  and  relevant  wavelengths. 

The  abscissa  in  the  plot  numbers  the  stations 
sequentially  starting  with  the  first  Duck  data  set  and 
terminating  with  the  last  Duck  data  set. 

As  is  seen  in  the  plots,  the  algorithm  in  some  cases 
does  a  good  job  of  giving  a  reasonable  value  of  the 
downwelling  k  for  certain  stations.  For  example,  for  the  490 
nm  k-value,  stations  1,  Z,  A,  and  6  are  fairly  well 
represented.  The  520  nm  k-value  is  closely  calculated  in 
stations  1,  2,  4,  and  6.  These  are  the  same  stations  as  for 
the  490  nm  result.  The  points  at  stations  3  and  5 
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show  far  more  erratic  behaviour. 


Station  3  here  is  the 


station  3  of  Duck  data  set  one,  and  station  5  is  the  second 
station  done  for  the  second  time  in  the  last  Duck  data  set 
Observing  the  irradiance  ratio  of  station  3,  in  Fig  12,  it 
appears  to  behave  in  an  unusual  manner  showing  a  large  'dip' 
in  these  ratios  as  a  function  of  depth 

No t hw  i  t hs t a n d i  n g  ,  the  Austin  and  Petzold  algorithm,  upon 
a  cursory  examination  of  the  calculated  and  measured  k’s, 
provides  a  seemingly  acceptable  way  to  find  the  k  of  a  water 
mass  from  orbit  or  from  an  airborne  platform.  However,  there 
appears  to  be  still  some  water  mass  features  that  may  inhibit 
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APPENDIX  F 
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